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ABSTRACT 



We have undertaken a new ground-based monitoring campaign to improve the 
estimates of the mass of the central black hole in NGC 4151. We measure the lag 
time of the broad H/? line response compared to the optical continuum at 5100 A 
^ ' and find a lag of 6.6lg § days. We combine our data with the recent reanalysis of 

■ UV emission lines by Metzroth et al. to calculate a weighted mean of the black 

^ ' hole mass, Mbh = {'^■^^-oij) ^ lO^M©. The absolute calibration of the black 
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hole mass is based on normalization of the AGN black hole mass - stellar velocity 
dispersion (Mbh — c*) relationship to that of quiescent galaxies by Onken et al. 
The scatter in the Mbh ~ cr* relationship suggests that rcvcrbcration-mapping 
based mass measurements are typically uncertain by a factor of 3-4. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 

1. INTRODUCTION 

Reverberation-mapping (Blandford & McKee 1982; Peterson 1993) is the process of 
measuring the time lag of variations in the emission lines of an active galactic nucleus (AGN) 
relative to variations in the continuum source due to light travel time effects. The time lag 
and the width of the emission line is directly related to the mass of the central black hole 
(Mbh)- Reverberation-mapping efforts over the past 15 years have led to the compilation of 
broad-line region (BLR) radius measurements and black hole mass estimates for 36 active 
galaxies (Peterson et al. 2004, 2005). Peterson et al. (2004) undertook a reanalysis of most 
of the reverberation-mapping data. The data were reanalyzed in a consistent way, paying 
particular attention to recent improvements in cross-correlation analysis, measuring the line 
widths in the variable part of the emission line, and accounting for the fact that changes in 
the mean luminosity of the source could cause emission-line time lags to vary over longer time 
scales than the reverberation time scale (Peterson et al. 2002). While the end result was to 
significantly decrease the amount of random and systematic errors in a very inhomogeneous 
database, there remained a few data sets for individual objects that were problematic in 
some way, making it difficult to determine either time lags and/or line width measurements 
with high confidence. NGC 4151, one of the brightest and best-studied AGN, is one such 
source. Previous monitoring studies of NGC 4151 include those of Maoz et al. (1991) and 
Kaspi et al. (1996). 

Our motivation for undertaking a new reverberation-mapping campaign on NGC 4151 
is twofold. Firstly, reverberation mapping is still the only direct method for estimating 
central masses in a broad population of active galaxies, and is more valuable still in that it 
does not depend on angular resolution. Reverberation results also provide the fundamental 
calibration for mass estimates based on radius-luminosity scaling relationships (Wandcl et al. 
1999; Vestergaard 2002, 2004; McLure & Jarvis 2002; KoUmeier et al. 2006; Vestergaard & 
Peterson 2006). Secondly, NGC 4151 is one of the best candidates for measuring the black 
hole mass by other means, and therefore affording a direct comparison between reverberation- 
based masses and other methods. Our specific goal in undertaking this new ground-based 
monitoring program is to significantly improve measurement of the lag, or mean response 
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time, between the continuum and the H/3 emission hne, the width of the H/3 hne in the 
variable part of the spectrum, and the black hole mass for NGC 4151. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Spectroscopy 

On every clear night between 2005 February 27 and 2005 April 10, we obtained spectra 
of the nucleus of NGC 4151 — a weakly-barred Sab Seyfert galaxy at z = 0.00332 — with the 
Boiler and Chivcns CCD Spectrograph (CCDS) on the MDM Observatory 1.3-m McGraw- 
Hill Telescope. We used a 5" slit with a position angle (PA) of 90°, and the typical exposure 
time was 1200 s. The spectra were reduced in the usual way using IRAF^ and XVista^. An 
extraction width of 17 pixels, corresponding to 12^75, was used in the spectral reduction. 
Each night had between two and four individual spectra of NGC 4151, which were averaged 
together to increase the signal-to- noise (S/N). 

Additional spectra were obtained with the Nasmith Spectrograph and Astro-550 580 x 
520 pixel CCD (Berezin et al. 1991) at the Crimean Astrophysical Observatory (CrAO) 2.6- 
m Shajn Telescope. Spectra were obtained through a 3" sht at PA = 90° and were reduced 
in the usual way. An extraction width of 16 pixels, corresponding to 11", was used in the 
reduction. 

Once the spectra were reduced, they were internally calibrated within each data set. 
Because the [O III] narrow hues are very bright in NGC 4151, we had saturation problems 
with the A5007 A line, and therefore scaled all the spectra using the A4959 A line. This 
was accomplished by employing the spectral scaling software of van Groningen & Wanders 
(1992). A mean spectrum derived from all of the data in each set was used as the reference 
spectrum. The individual spectra were compared to the reference spectra and scaled so as 
to minimize the residuals of the [O III] A4959 line in the difference spectrum (produced by 
subtracting the reference spectrum from the individual spectrum). This process results in 
a common [O III] A4959 flux of 3.76 x 10"^^ ergs s"^ cm^^. As the flux ratio of the [O III] 
A5007 to [O III] A4959 hne is 2.94, our measurement of the [O III] A4959 flux is consistent 



^IRAF is distributed by the National Optical Astronomical Observatory, which is operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under cooperative agreement with the NSF. 

^XVISTA was originally developed as Lick Observatory Vista and is now maintained in the public domain 
by former Lick graduate students as a service to the community. It is currently maintained by Jon Holtzman 
at New Mexico State University, and is available at http://ganymedc.nmsu.edu/holtz/xvista. 
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with the values of F([0 III]) = 1.575 x 10~^^ ergs s~^ cm~^ for the combined [O III] flux 
determined by Kaspi et al. (1996) and F([0 III] A5007) = 1.19 ± 0.06 x 10"^^ ergs s"^ cm'^ 
determined by Antonucci & Cohen (1983). 

The light curves were created by measuring the average continuum flux between the 
observed-frame wavelengths 5100-5130 A. To measure the H/? hne flux, we use observed- 
frame wavelength windows of 4770-4800 A and 5100-5130 A to estimate a linear continuum 
and then integrate the flux above the best fit continuum between 4810-4940 A. 

2.2. Photometry 

Photometric F-band data was obtained with the multicolor imaging photometer (MIP) 
at the 2.0-m Multicolor Active Galactic NUclei Monitoring (MAGNUM) telescope at the 
Haleakala Observatories in Hawaii (Kobayashi ct al. 1998a. b; Yoshii 2002; Yoshii ct al. 2003). 
The telescope was pointed first at NGC 4151 and then at two reference stars with {Aa,A6) 
= (16^1, 5'7) and (-11^3, -9^3), after which the telescope was dithered. The reference stars 
were previously determined to be non- variable (Minezaki et al. 2004). 

The images were reduced with IRAF in the usual way, with an additional small cor- 
rection applied for the nonlinearity of the detector. Aperture photometry with an aperture 
diameter of 8'.'3 was applied to the nucleus of NGC 4151 and the two reference stars, with 
sky flux measured between radii of 5'.'5-6'.'9. The nuclear flux of NGC 4151 relative to the 
reference star at {Aa,A5) = {I6'.l, 5'.7) was estimated to be the average of the relative 
fluxes for all of the dithering sets. The reference star flux was calibrated with respect to 
photometric standard stars taken from Landolt (1992) and Hunt et al. (1998). 

2.3. Intercalibration of Light Curves 

Each of these data sets involve different aperture geometries and seeing conditions, which 
result in different slit losses and varying amounts of starlight contamination. Therefore, 
they must be intercalibrated to flt each other before merging them into a flnal light curve. 
The MDM data set was the largest, and was therefore the anchor to which the CrAO and 
MAGNUM data were scaled. Following the discussion of Peterson et al. (1991), the CrAO 
and MAGNUM light curves were each compared to the MDM light curves to identify pairs 
of observations that were within one day of each other. A least-squares analysis was then 
performed on the sets of pairs of observations to identify the overall offsets, the results of the 
various monitoring geometries, between the MDM data set and the CrAO and MAGNUM 
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data sets. After removal of the offsets, the CrAO and MAGNUM hght curves were merged 
with the MDM hght curves. 

The full hght curves are given in Table 1. For our analysis, we binned all observations 
that occurred within a 0.5 day window. The statistical properties of the binned light curves 
are given in Table 2. Column (1) gives the spectral feature and column (2) gives the number 
of measurements in the light curve. Columns (3) and (4) are the sampling intervals between 
data points, measured as the mean and median, respectively. Column (5) is the mean flux 
and standard deviation of the light curve, not including a correction for host-galaxy starlight, 
which, following the methods of Bentz et al. (2006) for this slit geometry, is measured to be 
Fgai(5100A) = (2.11±0.18) X lO^"*^"^ ergs s""*^ cm" -2 A-^ Column (6) erives the mean fractional 
error, which is based on the comparison of observations that are closely spaced in time. The 
excess variance in column (7) is computed as 

F (1) 

where is the variance of the fluxes, 5^ is their mean-square uncertainty, and (/) is the 
mean of the observed fluxes. Finally, column (8) is the ratio of the maximum to the minimum 
flux for each lightcurve. 



3. DATA ANALYSIS 

3.1. Time Series Analysis 

To measure the lag between the continuum and the variable part of the H/3 line, we 
cross-correlate the H/3 emission-line light curve with the continuum light curve measured at 
5100 A. We use the Interpolation Correlation Function (ICF) method as described by White 
& Peterson (1994). The light curves and cross-correlation functions are shown in Fig. 1. 

To better quantify the uncertainties in the time delay measurement, we employ the 
model-independent Monte Carlo FR/RSS method, as described by Peterson et al. (1998), 
including modifications described by Peterson et al. (2004). For each single reahzation of the 
method, "random subset sampling" (RSS) is employed in which a light curve with data 
points is randomly sampled times without regard to the previous selection of each point. A 
data point that is selected M times will have its uncertainty reduced by a factor M'^^'^. "Flux 
randomization" (FR), in which a random Gaussian deviation based on the associated error 
bar, is then applied to each of the selected N points. This FR/RSS-altered subset of data 
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points is then cross-correlated as though it was real data. The peak of the cross-correlation 
function, Tmax, which occurs at a time lag Tpeak is determined, as is the centroid, Tcent; 
which is computed from those points near the peak with r > O.Srmax- A cross-correlation 
peak distribution (CCPD) for Tpcak and a cross-correlation centroid distribution (CCCD) for 
Tcent is built with a large number (A^ = 10000) of Monte Carlo realizations. We take the 
average value of the CCPD to be Tpeak and the average value of the CCCD to be Tcent, with 
uncertainties ATupper and ATiower such that 15.87% of the CCCD realizations have values 
T > Tcent + ^Tupper and 15.87% havc values t < Tcent ~ ■^'Hower, with similar uncertainties for 
Tpeak- Thcsc definitions of the uncertainty correspond to ilcr for a Gaussian distribution. 
We find an observed-frame lag time of Tcent = ^-^^-l/rl ^^J^ 'Tpeak — G.lOloieo days for 
NGC 4151. 



3.2. Line Width Measurement 

In order to calculate the black hole mass, we need to measure the width of the H/3 line 
in the variable (rms) part of the spectrum. The mean and rms spectra, shown in Fig. 2, 
were calculated using the full set of MDM spectra. To measure the line width, we first 
interpolate the rms continuum underneath the H/3 emission line by choosing continuum 
windows on either side of the line. The line width is then characterized by its full-width at 
half maximum (FWHM) and by the line dispersion (the second moment of the line profile) 
(Jiine, as described by Peterson et al. (2004). In order to characterize the uncertainties in 
each of these measurements, we use the following procedure, as described by Peterson et al. 
(2004). For a set of N spectra, we select N spectra at random, without regard to whether 
a spectrum has been previously selected or not. The N randomly selected spectra are then 
used to construct mean and rms spectra from which the line width measurements are made. 
This process is one Monte-Carlo realization, and a large number (TV = 10000) of these 
realizations gives a mean and standard deviation for each of the measurements of the line 
width. The 11/5 emission-line width measurements we determined from the rms spectrum for 
NGC 4151 are aune = 2680 ± 64 km s"^ and FWHM = 4711 ± 750 km s"^ in the restframe 
of NGC 4151. 



3.3. Black Hole Mass 

The mass of the black hole is given by 
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Mb„ = (2) 

where r is the emission-hne time delay, AV^ is the width of the emission hne, and G is the 
gravitational constant. The factor / depends on the geometry, kinematics, and inclination 
of the BLR and has been shown by Onken et al. (2004) to have an average value of (/) = 
5.5 when normalizing the AGN relationship between black hole mass and stellar velocity 
dispersion (Mbh — o"*) to that of quiescent galaxies. 

To estimate the mass of NGC 4151, we use Tccnt foi' the time delay and o"i;„c for the line 
width of H/?. We find Mbh = (S-l^ol) x 10"^ Mq. We would, however, like to remind the 
reader that all reverberation-based masses are likely to be uncertain by a factor of 3-4. 



4. DISCUSSION AND CONCLUSIONS 

It has been clear for many years that the BLR in AGN is stratified in ionization. Emis- 
sion lines arising from regions of more highly-ionized gas have shorter time lags, and lines 
arising from regions of less highly-ionized gas have longer time lags (Clavel et al. 1991; Pe- 
terson et al. 1991; Dietrich et al. 1993; Maoz et al. 1993). Additionally, it has been known 
that lines arising from more highly-ionized gas have broader profiles than lines arising from 
less highly-ionized gas (Osterbrock & Shuder 1982). If gravity dominates the kinematics of 
the BLR gas, then a virial relationship between the time lag and the emission line width. 
Ay oc T~^/^, will exist. In fact, this relationship has been found in NGC 5548 (Peterson & 
Wandel 1999) as well as a handful of other objects (Peterson et al. 2000; Onken & Peterson 
2002; KoUatschny 2003). 

We can investigate the presence of the expected virial relationship between time lag and 
line width in NGC 4151 by comparing our measurements with those of other monitoring 
programs for NGC 4151. A reanalysis of previous ultraviolet monitoring campaigns for 
NGC 4151 was recently undertaken by Metzroth, Onken, & Peterson (2006). Even though 
the International Ultraviolet Explorer {lUE) data set of Clavel et al. (1990) is undersampled, 
and the lUE data set of Ulrich & Horne (1996) is of a fairly short duration, the estimate of 
Mbh determined by Metzroth et al. (2006) of (4.14 ± 0.73) x lO'^M© is in good agreement 
with the value we obtain above. 

Figure 3 shows the lag and line-width measurements (given by Tccnt and ctuhc) for 
NGC 4151 from Metzroth et al. (2006) and from this work. We fit a power-law to the 
UV data from Metzroth et al. (2006) and the optical data from this work. Even with the 
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small range of values available for r and AV, we find that the best fit gives a slope of 
a = —0.76 ± 0.24, consistent with the slope of a = —0.5 expected if the virial relationship 
between the size of the BLR and the velocity of the gas exists. We do not include mea- 
surements from two earlier ground-based monitoring campaigns on NGC 4151, one in 1988 
(Maoz et al. 1991) and one in 1993-94 (Kaspi et al. 1996), for reasons we will briefly mention 
here, but are described in detail in Metzroth et al. (2006). The Kaspi et al. data set has 
a rather unconstrained lag time that is consistent with zero once the monotonic increase in 
the continuum and line flux is removed. The Maoz et al. data set, on the other hand, has 
a variable line spread function as well as low spectral resolution that causes the [O III] hnes 
to blend together. These factors make the hne width determinations suspect. 

Although the line widths in the Maoz et al. (1991) campaign appear to be unreliable, 
the time lag between the continuum and H/3 line is well deflned. The reanalysis of this 
data set in Peterson et al. (2004) gives a lag time of Tcent = 11-5 ± 3.7 days, which is 
quite a bit longer than the lag time found in this study. This discrepancy is not surprising 
given the variable nature of such objects. For comparison, the lag time of the H/3 line 
with respect to the continuum has been observed to range between 6.5 and 26.5 days in 
NGC 5548 as the ionizing flux level changes over time (Peterson et al. 2004). With such 
a difference in the measured lag times of the H/3 line in NGC 4151, we would expect that 
the flux would be 30% higher in the Maoz et al. campaign than in the current study, 
F{Hp) = (20.1 ± 2.5) X 10~^^ ergs s~^ cm^^ . In actuality, the flux during the Maoz et al. 
campaign, F{H(3) — (21.5 ± 3.4) x 10~^^ ergs s~^ cm~^, is 1.5(7 below the expected value. 
However, this is not outside the observed scatter for the H/3 lag times and luminosities 
measured for NGC 5548 (see Fig. 3 of Peterson et al. 2002). 

Table 3 hsts the black hole mass estimates for NGC 4151 calculated using a variety of 
data sets and methods. We show in bold face the weighted mean of the black hole mass 
estimated using the H/3 lag and line- width and the masses determined with the UV emission 
lines analyzed by Metzroth et al. (2006), a value of Mbh = 4.57+o;47 x IO^Mq, which we 
take to be the best reverberation-mapping mass estimate for NGC 4151. We can compare 
the mass determined using reverberation-mapping to the mass expected from the Mbh — ^* 
relationship. The stellar velocity dispersion of NGC 4151 was measured by Ferrarese et al. 
(2001) to be (7* = 93 ± 14 km s~^. Using the slope and zero-point determined by Tremaine 
et al. (2002), the Mbh — relationship predicts a black hole mass of Mbh = 6.2l|^xlO6M0, 
with both the scatter in the Mbh — c* relationship and the uncertainty in the measurement 
of (7* accounting for the large errors on the predicted black hole mass. This is about a factor 
of 7 smaller than the reverberation measurement. We caution, however, that the calibration 
of the reverberation-based black hole mass scale (Onken et al. 2004) is averaged over many 
objects. The scale factor / in eq. (1) depends on the structure, kinematics, and inclination 
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of the BLR (cf. Collin et al. 2006) and is expected to be different for each AGN. Use of 
the mean scaling factor (/) = 5.5 removes bias from the sample, in the sense that as many 
masses are underestimated as overestimated. Given that there is a factor of 3-4 scatter 
around the normalized AGN Mbh — cf* relationship, a single object deviating by a factor of 
seven is hardly inconsistent. 

In summary, the ground-based reverberation-mapping data presented here for NGC 4151 
supersedes previous optical data sets on this AGN. We have measured the lag of the H/3 line 
to the continuum fiux at 5100 A and found it to be G.G^Qg days. We have calculated the 
weighted mean of black hole masses determined using H/? as well as the UV emission lines 
analyzed by Mctzroth et al. (2006), which we find to be A.hltoTi x lO^M©. We also find 
that the data are consistent with a virial relationship between the size of the BLR and the 
velocity of BLR gas in NGC 4151. 
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Table 1. Continuum and H/3 Fluxes for NGC 4151 
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01 Ro 4- n 7n 

oi.DZ zt U. / U 


00/11 4- n /in 

zz. 41 It U.4U 




040 / .olZ 


08 1 /I -1- n Rf; 

Zo.l4 Zt U.DO 


OQ QQ 4- n OQ 

Zo.oo It u.zy 


IVi 


'^/I'^S 1/18 
040o.l4o 


OQ nn -1- n 07 
zy.uu zt u.z / 




TJ 

n 


Q/l QS 7i^n 
O4oo. / OU 


08 Qc: 4_ n RR 
Zo.oO ± U.DD 


OQ RQ 4- n Qn 

zo. oy it u.ou 


iVi 


70^ 
040y. ( ZO 


07 1 n -1- n RQ 
z ( .iu It u.oo 


OQ Q7 4- n OQ 

zo.o ( It u.zy 


iVi 


044U. ( 04 


OR nQ -1- n Rn 
ZO.UO It u.ou 


OQ 08 4- n OQ 

zo.zo It u.zy 


iVi 


Q/1/11 8Q8 
0441. oyo 


0/1 8Q -1- n c;8 

Z4.oy It U.Oo 


OQ 88 4- n Qn 
zo.oo It u.ou 


iVi 


OA An 7C/1 
o44Z. / 04 


OQ C^Q 1 n ^/l 
Ao.oo ± U.04 


OQ Q8 4- n Qn 
Zo.oo it U.oU 


iVi 


o/i/io son 

044Z.OZU 


0/1 1 c: -I- 1 1 
Z4.10 It 1.10 




tj 
n 


o44o. / 04 


00 Q8 4- n 1^0 
ZZ.Oo It u.oz 


OQ 0/1 4- n OQ 

Z0.Z4 It u.zy 


i\/r 

iVi 


0444. oyo 


00 Qn 4- n 1^1 
zz. yu It u.oi 


OQ n8 4- n AO 
zo. Uo It U.4Z 


>^ 


3445.391 


21.43 ± 0.47 


22.69 ± 0.41 


C 


3445.766 


21.49 ± 0.50 


22.36 ± 0.28 


M 


3446.332 


21.28 ± 0.47 


21.56 ± 0.39 


C 


3446.746 


21.96 ± 0.51 


21.85 ± 0.28 


M 


3447.828 


21.09 ± 0.69 




H 


3447.930 


21.66 ± 0.50 


20.83 ± 0.26 


M 


3450.785 


21.64 ± 0.50 


19.76 ± 0.25 


M 


3451.777 


21.62 ± 0.50 


19.34 ± 0.24 


M 


3451.988 


22.14 ± 0.16 




H 


3452.695 


21.65 ± 0.50 


18.99 ± 0.24 


M 
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Table 1 — Continued 





Fa (5100 A) 


H/5 A4681 


Observatory 


(-2,450,000) 


(10-15 ergs s"! cm-^ kr^) 


(10-1^ ergs s"! cm-^) 


Code'^ 


o40D. 100 


ZZ.U/ ± U.oi 


iy.Ui ± U.z4 


M 


o4oo.yiU 


zi.yu ± U.ol 


I T on _L n OO 

I I .o\) ± O.zo 


M 


o4oy. / lo 


Zi.ZO ± 0.49 


loin _i_ n OO 
lo.lO ± O.zd 


A /r 
M 


d4oU.oOo 


Zi. (4 ± U.OU 


1 o nn _i_ n oo 
io.UU ± U.zo 


A/T 
M 


d4oi. / 4o 


Zl.OI ± U.t)U 


I T o^ _i_ n OO 

I I .OO ± u.zo 


A/r 
iVl 


Q /I /^O 1 


OI on _i_ n /in 

zi.zy ± u.4y 


^ ^ p.^ _i_ n oo 
i / .Oi ± U.zz 


A/T 

M 


64:06. 6o I 


OO _i_ n txo 
zo.c)^ ± O.oz 


1 T _i_ n OO 
1 ^ . (4 ± U.dz 


/-^ 


o4d4.o/ 1 


OO no 1 An 
zo.Oz ± U.ol 


1 T CO _L n OO 

1 / .oz ± U.oz 


/-^ 


3464.809 


22.27 ± 0.52 


18.00 ± 0.23 


M 


3465.805 


21.93 ± 0.51 


17.52 ± 0.22 


M 


3466.793 


22.68 ± 0.52 


17.29 ± 0.22 


M 


3466.941 


21.87 ± 0.16 




H 


3467.895 


22.03 ± 0.51 


17.55 ± 0.22 


M 


3468.781 


20.66 ± 0.48 


17.70 ± 0.22 


M 


3469.398 


21.01 ± 0.46 


17.46 ± 0.31 


C 


3469.734 


20.80 ± 0.48 


17.56 ± 0.22 


M 


3469.941 


20.64 ± 0.11 




H 


3470.312 


20.80 ± 0.46 


17.28 ± 0.31 


C 


3470.785 


20.30 ± 0.47 


17.46 ± 0.22 


M 


3471.785 


20.66 ± 0.48 


16.90 ± 0.21 


M 



^The Julian Date listed is the midpoint of the observation 

^Observatory Codes: C = CrAO 2.6-m Shajn Telescope + Nasmith Spectro- 
graph; H = Haleakala Observatories 2.0-m MAGNUM Telescope + MIP; M = 
MDM Observatory 1.3-m McGraw-Hill Telescope + CCDS. 
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Table 2. Light Curve Statistics 



Sampling Mean 

Time Interval (days) Mean Fractional 

Series (T) T^edian Flux^ Error Fy^ R^^ 

(1) (2) (3) (4) (5) (6) (7) (8) 

5100 A 37 1.2 1.0 23.8 ±3.0 0.019 0.124 1.45 ±0.03 

H/3 34 1.3 1.0 20.1 ±2.5 0.012 0.121 1.41 ± 0.02 



Fluxes are in the same units as those used in Table 1. 



Table 3. Black Hole Mass Estimates 



Deseri])tion 




References 


Previous Ha + H/3 


1 qq+0.46 
-"-•<-> <J_0.46 


1 


UV Weighted Mean 


4 14+0.73 


2 


New R/3 


r nr. +0.89 
^■U'->-0.63 


3 


UV + up Weighted Mean 


-0.47 


3 


-^BH — cr^ Estimate 


62+°-^^ 

'-'•"^-0.30 


4,5 



References. — (1) Peterson et al. (2004); (2) Metzroth ct al. 
(2006); (3) this work; (4) Ferrarese et al. (2001); (5) Tremaine 
et al. (2002). 
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I I I I I I I I I L I I I I I I I I I I I I I I I I I I I I I L 

53440 53460 -20 -10 10 20 

JD (-2,400.000) Time delay (days) 

Fig. 1. — The left panel shows the lightcurves for the continuum region at 5100 A and for 
the H/3 A4861 line. Measurements taken within a 0.5 day bin have been averaged together. 
The continuum flux is in units of 10~^^ erg s~^ cm~^ and the H/3 flux is in units of 10~^^ 
erg cm~^ A~^. The right panel shows the result of cross-correlation with the continuum 
light curve; the top right panel is therefore the continuum autocorrelation function. 
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Fig. 2. — Mean and RMS of the MDM spectra in the observed frame of NGC 4151. The 
narrow components of the hnes have been subtracted from each spectrum before creating 
the RMS spectrum to avoid the residuals of the very strong [O III] hnes in NGC 4151. Some 
of the structure along the top of the H/3 broad hne is likely due to the imperfect subtraction 
of the narrow component of the line. The flux for both spectra is in units of 10~^^ erg s~^ 



cm 



A- 



-17- 



n 1 r 



Best fit slope: 
a= -0.76 ± 0.23 



• This work 

O Metzroth et al. 2006 



10 



cent 



(days) 



Fig. 3. — Rest-frame emission-line widths versus time lags for lines in NGC 4151. The filled 
circle is the data from this work. The open circles are UV emission lines from Metzroth 
et al. (2006). The solid line has a slope of — 0.76 ± 0.24 and is the best fit to the circles. The 
dotted line has a forced slope of —0.5, i.e. a virial relationship. 



